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Abstract 

The mean field like gauge invariant variational method formulated recently, is applied 
to a topologically massive QED in 3 dimensions. We find that the theory has a phase 
transition in the Chern Simons coefficient n. The phase transition is of the Berezinsky- 
Kosterlitz - Thouless type, and is triggered by the liberation of Polyakov monopoles, 
which for n > 8 are tightly bound into pairs. In our Hamiltonian approach this is seen as 
a similar behaviour of the magnetic vortices, which are present in the ground state wave 

functional of the compact theory. For n > 8, the low energy behavior of the theory is the 

* On leave of absence from ITEP, B.Cheremyshkinskaya 25, Moscow, 117259, Russia 
^permanent address 
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same as in the noncompact case. For n < 8 there are no propagating degrees of freedom 
on distance scales larger than the ultraviolet cutoff. The distinguishing property of the 
n < 8 phase, is that the magnetic flux symmetry is spontaneoously broken. 


2 



1. Introduction. 


In recent years topologically massive gauge theories, (TMGT) i.e. three-dimensional 
gauge theories with a Chern-Simons term [I] have attracted a lot of attention. 

In this paper we study the simplest model of this variety - the topologically massive 
17(1) gauge theory. The Lagrangian of this theory in the naive continuum limit is 

1 Tl 

L = -r—^F^F^ + —e^yx F^ A\ (1-1) 

where the charge g 2 has a dimension of mass and the Chern-Simons coefficient is a number. 
In a noncompact case this Lagrangian describes free massive particles with mass 2k = 
ng 2 /27t. We will be interested in the compact case. Our motivation to study this problem 
is twofold. 

First, compact three-dimensional QED without Chern-Simons term is confining due 
to Polyakov’s monopoles-instantons p. The coexistence of these monopoles with Chern- 
Simons term is an interesting question. It was studied to some extent in several papers 
p[ where it was argued that the monopoles are irrelevant in the presence of the Chern- 
Simons term, the theory looses its confining properties and behaves in all aspects as a 
non-compact theory. On the other hands these arguments are not very rigorous and the 
question of the relevance of the monoples is not completely settled. 

The second issue arises from the connection between the 2 + 1 topologically massive 
gauge theory with compact 17(1) group and an induced 2 dimensional XY model on the 
space - time boundary of the 2+1 dimensional manifold. This model arises because the 
Chern-Simons term in the TMGT Lagrangian ( |1.1|) is not gauge invariant under the gauge 
transformations which do not vanish on the boundary. These gauge degrees of freedom 
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become physical and have nontrivial dynamics on the boundary. One can show that the 


induced two-dimensional action has the form 



J d 2 x di4>di4> 


( 1 . 2 ) 


This is the action for the massless two-dimensional scalar field (f) which is a parameter 
of a gauge transformation in the original three dimensional theory. The question of 
whether the three dimensional theory is compact or not is of crucial importance. In 
the noncompact case the gauge parameter varies from — oo to oo and the theory on the 
boundary is the theory of free massless bosons - the conformal field theory with c = 1. 
In the compact case with gauge group 17(1), the gauge parameter lives on a circle S 1 
and eq. (|1.2j ) defines the XY model which has the famous Berezinsky-Kosterlitz-Thouless 
(BKT) phase transition [j|]. The difference between the XY model and the free scalar field 
is the existence of singular field configurations (vortices). The XY model has two phases. 
At large n the vortices are bound into dipolar molecules and do not affect long distance 
physics. This phase is described by the conformal c = 1 theory. At small n the vortices 
are in the plasma phase. They dominate the statistical sum and lead to emergence of 
finite correlation length. The critical value of n for the action (|1.2|) is n = 8. Thus, due to 
the connection between the XY model and the compact U( 1) TMGT one is lead to think 
that the three-dimensional theory should also have two phases. This phase transition and 
existence of a new phase in the topologically massive gauge theory as well as some related 
issues, for example n —> 1 /n duality in both theories (R —> 1/R in a standard notation 
in string theory) have been discussed some time ago by one of the authors ||. However 
the complete picture of the phase transition at small n was not clear at that time. 

In this paper we study the topologically massive compact 17(1) gauge theory using 
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the gauge invariant variational approach developed in In ref. |7j this method 

was used to study compact QED 3 without the Chern Simons term. These results are in 
agreement with the well known scenario of confinement due to magnetic monopoles. In 
the framework of the hamiltonian approach of 0 the monopoles appear as topologically 
nontrivial configurations of the gauge group when the trial vacuum wave functional is 
projected onto the gauge invariant subspace of the Hilbert space. The advantage of 
this approach is a straightforward generalization in the case of non-zero Chern-Simons 
term. As opposed to that, in the path integral approach in the presence of the non-zero 
Chern-Simons term, the monopoles are not solutions of the classical equations of motion 
anymore. One has to consider complex-valued solutions of the equations of motion and 
their meaning is not completely clear [3|. 

We show that the model has two phases depending on the value of the Chern-Simons 
coefficient n. The existence of these two phases is related to the different behaviour of the 
monopoles. At n > 8 the monopoles are irrelevant for the infrared physics and the theory 
behaves as in the noncompact case and induces conformal c = 1 model on the boundary. 
This is also in a qualitative agreement with the results of ref. |3|. 

For n < 8 the monopoles condense in a new ground state. In this phase they are 
important for the infrared physics and the induced theory on the boundary is a deformed 
c = 1 model corresponding to the plasma phase. Let us note that even though there is a 
monopole condensate in this phase, there is no confinement for any nonzero n. 

The transition between these two phases indeed takes place precisely at n = 8, which is 
a strong argument in support of the idea that all phenomena in induced two-dimensional 
dynamics have they counterparts in an original three-dimensional theory. 

The organisation of the paper is the following. In the next section we consider the 
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Hamiltonian formalism for the topologically massive gauge theory. First we discuss the 
noncompact case and then modify it for the case of compact 17(1). We introduce the 
gauge invariant trial wave functional. In section 3 we discuss some properties of this wave 
functional and the interpretation of the calculations in terms of gases of magnetic and 
electric vortices. In section 4 we calculate the average energy in this trial functional and 
minimize it with respect to our variational parameters. We will demonstrate that for large 
n the long distance properties of the ground state are the same as in the noncompact case 
and then will investigate the case of small n. In conclusion we discuss the properties of 
the new n < 8 phase. 

2 The model and the setup. 

2.1 The noncompact limit. 

Let us start with setting up the Hamiltonian formalism for the theory. To do that consider 
first a noncompact theory. It is described by the following Hamiltonian 

ff=I[£? + B 2 ] (2.3) 

augmented by the Gauss’ law constraint 

C(x) = di(Ei - 2 KdjAj) = 0 (2.4) 

The commutation relations satisfied by the fields are 

[Ai(x), Aj(y )] = 0; [E^x), Aj(y )] = iS i:j S(x - y ); [E^x), Ej(y )] = 2iKe ij 8{x - y) (2.5) 

The electric held E^ can be represented in terms of the momentum, canonically conjugate 
to the vector potential A t as 

Ei = —n* + K6ijAj (2-6) 
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This Hamiltonian formulation follows directly from the Lagrangian (1.1). 


The Gauss’ law operator C(x) generates time independent gauge transformation. The 
elements of the gauge group are the operators 


Utp = exp 


i j d 2 x(p(x)C(x)^j = exp ji J d 2 xdi(j){Jli{x) + K€ijAj(x))^j (2.7) 


ft is straightforward to check, that both E t and B commute with the Gauss’ law, and 
are therefore gauge invariant operators. The physical Hilbert space of the theory contains 
only states which are invariant under the action of U 


|T >= T > (2.8) 

It is most convenient for our purposes to work in the field basis. Then, defining 


$[ 4 ] =< 41 $ > 


(2.9) 


we find 

< >= exp ji/c J d 2 xdi(fr(x)eijAj(x ) j <f>[4 + <9,0] (2.10) 

A gauge invariant state can be then constructed from an arbitrary state $ by averaging 
over the gauge group 

$[4] = j D(j) exp j d 2 xd i (j)(x)6ijA j (x)^ $[Ai + di(j)} (2.11) 

In the case of noncompact theory the integral over 0 can be performed, and one can give 
a more compact representation of a gauge invariant state as 

rv 

exp{m J -£AiB}$[B] (2.12) 


However, we are going to deal with the compact theory, in which case the representation 
eq.([2.11|) is more helpful. 
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2.2 Compactifying the gauge group. The vortex operator. 


The next step in our discussion is to compactify the gauge group. As discussed in some 
detail in |7| this means, that we have to enlarge the gauge group so that it includes 
operators which create pointlike magnetic vorices of integer strength. These operators 
satisfy the following commutation relation 

0<7r If 

Vi(x)B{y)V k (x) = B(y) +-A 2 (cr - y) (2.13) 

9 

with g - a dimensional constant, determining the radius of compactness and k - arbitrary 
integer. Naively this just means including in the gauge group the gauge transformations 
of the form eq. (|T7|), corresponding to singular functions 

Tl 

4>(x) = — 9(x — xq) (2-14) 

9 

where 6 - is the planar angle, and with the understanding that the derivative of 0 in eq. 
( |2.7|) should be taken modulo 27T. For future convenience let us introduce a special symbol 
to denote this derivative 

A i9{x) = di[9 (x)\ mod2n = (2.15) 

x A 

That is, these derivatives do not feel quantized discontinuities in (p(x). 

In other words, to compactify the theory we must limit physical Hilbert space to states 
which are eigenstates of the operators 

V(x) = exp | ^ J d 2 yAi9(x - y)[Ei(y) - 2K6 ij A j }\j (2.16) 

with a unit eigenvalue. For k = 0 this procedure is straightforward and was described in 
[f]. However, for a nonzero k one should be a little more careful. The point is, that the 
operators defined in eq. (|2.16|) do not commute with noncompact gauge transformations, 






and also do not commute between themselves, but rather satisfy 


V(x)V(y) = V(y)V(x)e iKa ^'> 


(2.17) 


where a(x — y) is a c-number. Obviously, if t%a 2nm, one can not impose the condition 


y(x)|^ >= |t > simultaneously for all points x. Fortunately, eq. ( |2.13| ) does not define 
the operator V completely, and we can take advantage of this freedom. Let us consider 
the following operator 


V c (x) = exp j^y d 2 y[Aid(x - y)Efiy ) - 2ne ij d i 6{x - y)A j \'j (2.18) 

Note, that the function multiplying Aj in the exponential is now a bona fide derivative. 
Since the planar angle is defined relative to some direction, the function 6 has a discon¬ 
tinuity of 2n along some curve C, which starts at the origin and goes to infinity. Its 
derivative therefore has the form 


diO(x ) = A id{x) + 2'Keijh ( j (x)8{x — C) (2-19) 

where nf(x) is the unit tangent vector to the curve C at the point x. The operator V(x) 
defined in eq. ( |2.18| ) still satisfies eq. ( p,13| ), but now is invariant under the noncompact 
part of the gauge groupQ, and satisfies 

VoAnVcM = VcMVcX (2.20) 

where is the number of intersections between the curves C\ and C^- It is clear 

therefore, that if we choose the radius of compactness so that 

47 tk = ng 2 (2-21) 

*To see this, note that the A - dependent term in the exponential, integrating by parts can be rewritten 
as f (f>B 1 which is explicitly gauge invariant. In principle one could worry about the surface term, which 
we omitted / x 2 _ too ds l ipCijAj. In the present case however, it is harmless for two reasons. First, because 
the vector potential itself is massive, and therefore vanishes at infinity; and second, since for all practical 
purposes it is enough to consider transformations with nonzero number of vortices and antivortices, but 
with net zero vorticity, and for those the function <f> itself vanishes at infinity. 
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with arbitrary integer n, the operators Vc(x) commute at all points and for arbitrary 
choice of the set of curves C. Only in this case the compact theory will not depend on the 
choice of the curves, and will therefore be rotationally and translationally invariant. We 
see, therefore that for nonzero n the radius of campactness must be quantized according to 


eq. ( 2.21 ). This is, of course in complete agreement with the fact that in nonabelian Yang 
Mills theories, where the gauge group is necessarily compact, the ratio of the coupling 
constant and the Chern - Simons coefficient is quantized in the same way. 


The expression for the vortex operator V is simpler on a subspace of the Hilbert space 
invariant under the noncompact part of the gauge group. On these states the Gauss’ law 
is implemented exactly. Integrating the second term in the exponential in eq. (|2.18|) by 
parts and using the Gauss’ law, we get 


V s = exp{i J Si(y)Ei(y)} (2.22) 

with Si(x) = 2tt / geijfi^(x)5(x — C ). We find it more convenient to label the operator V 
by the vector function ,s t rather than by the curve C. They are of course in one to one 
correspondence. 


There is one subtlety related to the algebra of the operators V s for odd n. Although 
the vortex operators commute for any n, their group multiplication properties are not 
completely trivial. 

V S1 V S2 = e^ nN ^V sl+S2 (2.23) 

where N(l,2) is again the intersection number of the two curves. For even n the extra 
phase factor is always an integer of 27 t. In this case one can construct physical states by 
requiring V s x l ; = T for all s. For odd n the phase factor can be an odd integer of 7r. We 
therefore can not require that the wave functional be invariant under the action of all V s . 
This however can be remedied in the following way. Let us define for future convenience 
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the densities p and er by 


9 rj S' Q 

—eijdiSj = p, —OiSi = a 
Ztt Ztt 


(2.24) 


With our definition of s*, both p and a satisfy 


/ d 2 xp(x ) = 27rni, / d 2 xa(x) = 2 nn 2 
J S J s 


with integer ri[ and 5 for any finite area S'. Then the field r 


(2.25) 


d 2 r = 2i rcr 


(2.26) 


takes integer values. 


Now, note that the phase factor in eq. (|2.23|) can be represented as 
exp{*7rniV(l, 2)} = exp {in J 6ijs\s 2 } = expjy—— J [r 1 p 2 — pV 2 ]} 


Due to the property eq. (|2.25|) this can be further rewritten as 


ex P{^ J t r V 2 + P 1 ^ 1 }} = exp{i^ / [(r 1 + r 2 )(p 1 + p 2 )]} 


It is clear now, that a modified vortex operator, defined by 


= exp{i J SiEi + J r(s)p(s)} 


has a simple group multiplication rule 


(2.27) 


(2.28) 


(2.29) 


Ksi+S2 


(2.30) 


One therefore can require consistently, that these operators for all s t leave physical states 
invariant. In the following we will, however disregard this subtlety and for simplicity use 
the vortex operators defined in eq. ( |2.22| ). Thus our calculations in the way they are 
presented below are valid only for even n. 
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Before we continue, we wish to remark that the densities p and a as defined in eq. 
( 2.25|) have a very clear physical meaning. The vortex operator V s has the following 


commutation relations with electric and magnetic fields 

O tt 

V(x)B{y)V s {x) = B(y) -f p(y) (2.31) 

9 

(x)d i E i (y)V s (x) = diEi(y) + ngp(y) 

V} (x) €ij diEj (y)V s (x) = e i:j diEj (y) + nga (y) 

Therefore the operator V s when acting on any state creates magnetic vortices (with mag¬ 
netic flux quantized in units of 2 n/g) with the density p, and also creates electric vortices 
(with vorticity quantized in units of ng) with the density cr. 

2.3 The Hamiltonian. 

As just noted, the operators B and E t do not commute with V. The Hamiltonian of the 
noncompact theory should, however be invariant under the complete compact group. It 
must therefore be slightly modified from the simple form of eq. ( |2.3|) . The most natural 
way to do it, is to use gauge invariant operators which in the naive continuum limit reduce 
to the standard B 2 and E 2 terms in the continuum hamiltonian. For the magnetic part 
of the energy density we take therefore 

H b (x) = -?A[1 - Ree i9ma2B ^} (2.32) 

g-m 1 

Here a is the ultraviolet regulator which has a dimension of distance, and has the meaning 

of the lattice spacing.0 The constant m is an arbitrary integer, m < n. 

Hn general, in the following wherever ultraviolet regularization is needed we assume lattice regular¬ 
ization with lattice spacing a. Accordingly, all derivatives should be understood as symmetric lattice 
derivatives etc. 
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As opposed to the case k — 0, the electric part of the Hamiltonian should also be 
modified. Again, guided by the naive continuum limit and gauge invariance, we write for 
the electric part 

H e (x) = 9 1 ° Y [1 - Ree^ x)Ei{x) } (2.33) 

2?r i=l,2 

The vector function ai(x) is defined on the links of the lattice in terms of the unit vector 
fii(x) parallel to the given link as 




2n a 
9 


eijhjix) 


(2.34) 


In order for He to commute with the vortex operator V, the integer l must be a divisor 
of n: n/l = integer. Note, that for / = 1, the operator He(x ) is equal to the sum of two 
vortex operators. In that case, it is trivial on all physical states. The same is true for Hb 
for m = n. In that case, using the Gauss’ law one can see that Hb is given entirely in 
terms of a vortex operator V s where the function s* corresponds to the contour C which 
is the boundary of the elementary plaquette on the lattice. Therefore it is obvious, that 
in order for the dynamics of the theory defined with this Hamiltonian to be the most 
nontrivial, one should choose the largest possible l and the smallest possible m. However, 
it turns out, that to avoid subtleties similar to the ones encountered for odd n, it is 
convenient to choose both m and n/l to be even numbers. We will therefore take m — 2 
and l = n/2 in the rest of this paper. 
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3 The Variational wave functional. 


Now it is straightforward to write down a general wave functional, which is invariant 
under the whole compact gauge groupjj] 

^[A] = J D.SiD(j) exp jm J d 2 x ^(s^x) - d i <p(x))A j (x) + e ij d i (j)(x)s j (x)\ j <f> [A-<%</>-A 

(3.35) 

Our purpose is to find a wave functional of the vacuum of this theory using a gauge 
invariant generalization of a gaussian variational approximation. We shall minimize the 
VEV of the energy on the following set of explicitly gauge invariant states: 

^[A] = J DsiDcfiex p ^[sieijAj - d^e^Aj + e^d^Sj] - ^Af^G^Af’ 3 ^ (3.36) 

For convenience we have introduced the following notation 

Af (x)Mx) = Ai(x) - di(j>{x ) - Si(x) (3.37) 

and have switched to the matrix notations, so that 

diMbi — J d 2 xd 2 yai(x)Ad(x — y)bi(y) (3.38) 

In the next section we are going to calculate the expectation value of the Hamiltonian 
on this set of states, and functionally minimize it with respect to G(x — y), which is our 
variational parameter. 

3.1 The noncompact limit. 

We start with solving the noncompact limit of the theory. This serves to illustrate the 

method, and also to get a clear idea of what to expect of the variational function G(x). 

^The last, c-number, term in the phase factor appears, since the two components of the electric field 
do not commute, and is readily obtained using Baker-Campbell-Hausdorff formula. 
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Our variational state now is given by eq. (|3.36|) , but without the integration over s,;. The 
Hamiltonian is the simple quadratic Hamiltonian of eq.(2.3). The exact ground state 
wavefunctional in this case is therefore gaussian. So our variational calculation in this 
case will give the exact result. The expectation value of any gauge invariant operator is 
given by the following path integral 

< 0[A\ >= Z — j D(f)DAi exp j ~AfG _1 Af + wdifytijAj j 0[A\ exp j 

(3.39) 

The integrals over Ai and (j) are both Gaussian, and are easily performed. Calculating the 
relevant averages we obtain 


< E 2 (x) >= 


< B 2 (x ) >= 
1 f d 2 p 


d 2 p 


■p 2 [G 1 (p)+n 2 G(p)Y 


2 J (2tt) 2 

G _1 (p) + k 2 G(p) + Ak 2 [G^ 1 (p) + k 2 G(p)}~ 1 


2 J (2t r) 2 

Here G(k) is the Fourier transform of the variational function G(x). 
The expectation value of the energy is 


(3.40) 

(3.41) 


2V~ l <H>=- 


d 2 p 


G (p) + k G(p) + (p + 4 k )[G (p) + k G(p)] 


mi 


2 J (2t r) 2 L 

Minimizing this expression with respect to G(p) we obtain 


(3.42) 


G 1 (p) + n 2 G(p) = \Jp 2 + 4k 2 


(3.43) 


As is clear from this calculation, the function 


D{p) = [G {jp) + k G(p)} 


-i 


(3.44) 


plays the role of the propagator. For example, the equal time propagator of magnetic 
held is given by 


< 


B(x)B(y) >=\J A£-e‘^-’\2n)Y[G-\p) + K 2 G(p)]~' (3.45) 
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As expected, D(p) in the noncompact case comes out as the propagator of a free massive 
particle with mass m 2 = 4k 2 . 


3.2 The norm of the state in the compact theory. 

As a preamble to calculating the energy expectation value, let us find the norm of the 
state eq. (|3.36| ). 

Z = J DsiDcpDAi exp + Af s G~ l Af s ] - i/c[(s* - d^e^Aj + e^diSj] j 

(3.46) 

The integrals over A t and over the noncompact part of the gauge group, </> are gaussian. 
After performing them we are left with the integral over the descrete variable s;. 


Z — ZaZ^Zg 


(3.47) 


with 


Z a = det (7T G) 


(3.48) 


Z# = j D(f) exp {--/ di(j)D 1 d i 4> j = det 


An 


D 


-l 


- 1/2 


z s = f D.Si exp {-\ f d P 


Ai (2vr) 2 
where D(p) is dehned in eq. (|3.44|). 


Si(p)[$ij ~ r]D 1 (p)s j (-p) + n 2 s i (p)^AD(p)s j (-p) 


p 


p* 


A more convenient and physically intuitive representation of the integral over Sj can 


be given in terms of the variables p and a dehned in eq. (|2.24j) . 


Z s — Z p Z a 


(3.49) 


with 


^p = 


jDpexp | ~ j -A£-p(p)± D -\p)p(-p) 


g 2 J (2vr) 2 p 


(3.50) 
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(3.51) 


Z a = j Da ex p \ - 7 tuk / ^ a(p) .D (p)a (--p) 


(2tt) 


P 


One should keep in mind, that the measures Dp and Da are not the same as the 


measure for the functional integration over a free held. Rather, due to eq. (|2.25| ) both 
these path integrals are equivalent to statistical sums for classical interacting gases: 


oo n+ n- 


z P = 


exp 


y nn(« 2 ) n++n ~dx a dxp 

n_|_,n_= 0 a=l /3=1 


(3.52) 


(Y u { x a - x a') + Y U ( X P ~ X P') ~Y U ( Xa ~ X d)\ + ^ n+ ^ n 


P,P' 


a,(3 


oo n+ n- 

Z a = Y II Yl(a~ 2 ) n++n -dx a dxp 

n_|_,n—=0 a=l (3=1 


(3.53) 


{ ^ ( n + + n _) 

-\Y V \ X <x ~ X <*') + Y V \ X P ~ X P') - Y V \ X <x ~ X p)\ H- n — v (°) 

a,a 1 /3,/3' a,13 ~ 

The interparticle interaction potentials u and v are given by 


2 vr f d 2 p 1 x , 
u(x) = — 5 - / ———L> (p)cos(px) 


r 


v{x) = 27rm J 


(27t) 2 p 2 
d 2 p 1 


(3.54) 


D(p) cos (px) 


(27r) 2 p 2 ' 

Since u{ 0) and u(0) are singular, the last terms in the exponential in the equations 
Q3.53| . |3.54] ) should be understood, as usual in the regularized sense, that is at finite UV 
cutoff: w( 0 ) (u( 0 )) should be substituted by u{x = a) {v{x = a)). 

Remembering the interpretation of p and a, discussed in the previous section, we shall 
refer to the statistical mechanical systems defined by eqs.( [3.53|) and Q3.54D as gases of 
magnetic vortices and electric vortices respectively. 

The calculation of any expectation value (such as the expectation value of the energy) 
in our trial wave functional should now proceed in the standard way. The integration over 
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the vector potential A* and the noncompact part of the gauge group is always Gaussian 
and therefore trivial. After performing this integral we will be invariably faced with the 


problem of calculating certain correlators in the two vortex ensembles, eq. ( 3.53 , 3.54 ). It 
is therefore worthwhile to try and understand some general features of these ensembles. 


First note, that the cases k = 0 and k ^ 0 are essentially different. For k = 0 (or 
equivalently n — 0), the electric vortex partition function becomes completely trivial - 
the interaction potential vanishes. This is of course to be expected, since we know that 
the electric vortices do not play any role in the theory without the Chern - Simons term 


Another crucial difference is that the behaviour of the magnetic vortex interaction 
at large distances is significantly changed. At n = 0 the function = G~ l (p) for 

the best variational state vanishes at zero momentum |7|. Consequently, the interaction 
potential between the magnetic vortices is short range. For k ^ 0 on the other hand this 
cannot happen, since the expression D~ 1 (p) = G~ 1 (p) + k 2 G{j>) is bounded from below by 
2 k. Therefore the interaction between the vortices is logarithmic at large distances, and 
whatever G is, we are dealing with the Coulomb gas. Obviously, for large k the gas will 
be in the molecular phase, and one expects that it will have no effect on the large distance 
physics. This is in agreement with the general arguments of ref.||, that monopoles should 
be irrelevant in a Chern - Simons theory. However at smaller k the Coulomb gas will be 
in the plasma phase, and will certainly affect physics. The effective temperature of the 
Coulomb gas of magnetic vortices in our case is T m = 2g 2 /2 k = Anjn, and one therefore 
expects this change in behaviour to set in at n — 8. 

For nonzero n the electric vortices behave in a similar way. Taking the noncompact 
expression D~ x (p = 0) = 2k, we find that at large distances the interaction between the 
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electric vortices is also logarithmic, and in fact has the same asymptotics as the magnetic 
vortex potential. The effective temperature is also T e = Air jn and one again expects the 
Berezinsky-Kosterlitz-Thouless like phase transition at n — 8. 


3.3 The magnetic vortex gas. 


In the following we will need to calculate correlation functions of the vortex densities. 
To facilitate this we use the standard trick [:2|], j8| to rewrite the partition function of a 
classical gas in terms of a path integral over a scalar field. First consider the magnetic 
vortices. The exponential factor in eq. ( |3.53 ) can be rewritten as 

-.2 


i 2(n++n - } J D X exp{-^ X P 2 Dx + ipx} 


(3.55) 


The summation over the number of vortices (or integral over p) can be performed exactly 


Dpexp{ip(x)x(x)} = <5(exp{ix(a;)} - 1) 


(3.56) 


This gives 


Z n = 


f P 2 

J D\ .:xp!- — X:) rD 4 ;}n.,5!cxp!)\f.r)} - 1} 


(3.57) 


= / D X Da exp 


47T 2 


Xp 2 Dx + i J 2a(x)(cosx(x) - 1) 


In general the integration over the lagrange multiplier field a is nontrivial. However, 
when the density of vortices is small, one can perform the integration over p in the dilute 
gas approximation, that is sum only over those configurations which have only at most 
one vortex or antivortex in every point in space. This is equivalent to substitute for a a 
constant 

Z p = J Dxexpl—^xP 2 Dx + 2a- 2 cos X (x)\ (3.58) 
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To calculate the correlator of p in the dilute gas approximation one can add ipJ to the 
vortex free energy, and calculate functional derivatives of the resulting partition function 
with respect to J at zero J. A simple derivation gives 

< p{x)p(y) >= u^(x - y)- < u~ l x{x)u~ l x{y) > (3.59) 


The propagator of x is easily calculated. At weak coupling the interaction in eq. (|3.58|) 
is very small (< cosy ><< 1). To first order in the interaction the only contribution 
to the propagator comes from the tadpole diagrams. This is easily seen by rewriting the 
cosine potential in equation fl3.58|) in the normal ordered form 


cosy =< cosy >: cosy := za : cosy : 


(3.60) 


One should be a little careful in the definition of the normal ordering. Ordinarily the nor¬ 
mal odering would be performed relative to the free theory with the propagator u(x — y ). 
In the present case however the free propagator at small momentum behaves like a prop¬ 
agator of a massless particle u(k) —■>k->o k 2 , and the bubble integral which enters the 
calculation of < cosy > is infrared divergent. This problem can be overcome by perform¬ 
ing normal ordering relative to a massive theory. This can be done selfconsistently, by 
including the quadratic term in : cos(y) : into the free propagator. In this approximation 
the propagator of y is 


d 2 xe lkx < X (x)x(0) >= — 

1 i x 


u~ 1 (k) + 2z 


(3.61) 


with z determined selfconsistently by 


z = a 2 < cos(x) >= a 2 exp{-- < y 2 >} = 


(3.62) 


a 2 exp 


1 r d 2 p 

2J (2^p 


■^p 2 D(p) +2 z 


i -i 
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The correlator of the vortex densities is then 


KM = jdv* < P (x ) P (o) >= t+ 2 : u(k) 


(3.63) 


The existance of the critical point n — 8 is straightforward to see in this approximation. 
Using the result for the noncompact theory -D(O) = 1/2 k, and anticipating the fact that 
the infrared assymptotics of the propagator is the same in the noncompact theory, we can 
rewrite eq. (|3.63|) as 

z = a~ 2 (za 2 e~ c ) n/8 e~ flp (3.64) 


Here the constant c is the chemical potential of the pure Coulomb gas. It’s exact value is 
not important, but one has to remember that it is of order one and positive. The chemical 
potential /i p measures how different the magnetic vortex gas is from the Coulomb gas and 
is defined by 


dp 


d 2 p 

W 2 


2vr 2 


p 2 D(p ) + 2z\ - p 2 + 2 zj 


-l 


(3.65) 


The chemical potential p p depends on z very weakly for small z. Also, at weak coupling 
the chemical potential is very big. Taking, for orientation the noncompact result for 
D(p), we find that p p oc Under these circumstances it is easily seen that for n < 8 


the solution for eq. (|3.65|) exists, and is indeed parametrically much smaller than the 
ultraviolet cutoff a~ 2 . 


z oc a 2 e s -n^p 


(3.66) 


For n > 8 one can check that the selfconsistency equation eq. ( |3.63|) has only the 
trivial solution z = OPl. 


§ Even though eq. ( 3.66 ) for n > 8 seems to give z oc 0(a 2 ) this is not correct, since this equation 
was derived under assumption that z is small. 
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Importantly, the use of the noncompact expression for D(p) is not at all crucial for 
this derivation. Remember, that quite generally D(p) < 1/2 k. The chemical potential 
for the magnetic vortices is therefore always larger than the chemical potential in the 
Coulomb gas that corresponds to the XY model. The dilute gas approximation is known 
to work reasonably well for the XY model ||], and we are assured therefore that the 
dilute vortex gas approximation for the magnetic vortex gas should be reliable for any 
variational function G(p). 

3.4 The electric vortex gas. 

While the magnetic vortex gas can be consistently treated in the dilute gas approxi¬ 
mation, the behaviour of the electric vortex gas is more complicated. A dilute electric 
gas approximation is easily formulated along the same lines. This results in a partition 
function: 



(3.67) 


The “bubble summation” then gives the propagator of the field 0 



(3.68) 


with ( determined selfconsistently by 



(3.69) 


This can be recast into the form 



(3.70) 


with the chemical potential fi a 
1 f d 2 p 

^ = 2 J (2YP _ 


/ 


2'Ktnn 


1 



-il 


(3.71) 
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One should however ask oneself, under what circumstances is the dilute gas approx¬ 
imation valid. A necessary condition for this is that the interaction energy between the 
vortices is not negligibly small even at very short distances of order of the ultraviolet 
cutoff. It is clear therefore, that the diluteness of the electric vortex gas holds only for 
some restricted set of the variational functions D. Let us take, for example the noncom¬ 
pact form for D(p ), so that for small momenta D(p) = l/2tz and in the ultraviolet region 
D(p) oc (p 2 ) -1 / 2 . In this case the interaction potential between the vortices is logarithmic 
at large distances (x > 1/iz) but at short distances is very weak v(x ) oc ntzx. The vortex 
gas will therefore be dilute on all distance scales only if n > (iza)^ 1 , that is either for very 
large values of n, or in a theory where the mass of the photon is close to the ultraviolet 
cutoff. The former case corresponds to extremely weak coupling, while the latter to the 
limit of pure Chern - Simons theory without the Maxwell term. On the other hand for 
finite n and na —> 0 (naive continuum limit), the electric vortices therefore hardly feel 
the presence of each other if the distance between them is smaller than 1 / tz. At these 
distance scales the gas will therefore not be dilute at all. To wit, in this case the chemical 


potential ji c defined by eq. (|3. 71|) is large and negative. 


Another situation in which the dilute gas approximation is valid, is for those variational 
functions D(jp ) which have a constant value D[jp) = 1/2 tz for all momenta larger than some 
scale /i, such that fia « 1. In this case /i CT is close to zero, and one is back to the case 
of pure Coulomb gas. For these variational functions the behaviour of the electric vortex 
gas is very similar to the behaviour of the magnetic vortex gas. Both have BKT phase 
transition at n — 8, both are in the molecular phase for n > 8 and in the plasma phase 
at n < 8. 


The question is however, is there any reason to expect, that the best variational 
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propagator will behave in this fashion. Surprising as it may seem at the first glance, 
the answer to this is positive. In fact, a little thought convinces one that it is almost 
unavoidable in a compact theory with finite radius of compactness 1/g. The point is the 
following. As discussed in Section 2, imposing compact gauge invariance, among other 
things has an effect of imposing the following condition on the magnetic field 

>= |T > (3.72) 

This means, that a gauge invariant state T > has nonzero projection only on states with 
quantized eigenvalues of the magnetic field B: B(x) = We know, however that in 

the noncompact vacuum < B 2 >oc a~ 3 . This means that the natural scale for the magnetic 
held in the vacuum state of the noncompact theory is B oc a~ 3//2 . Therefore if one projects 
a state with the ultraviolet properties of the noncompact vacuum onto a gauge invariant 
state, only the contributions of B = 0 states will survive and the magnetic energy will 
vanish. Since the expectation value of the magnetic field in a noncompact state is just 
given by < B 2 >= / |-§ p 2 D(p), a nontrivial compact dynamics can be described only by 
states whith a much larger value of the “propagator” D(p). In fact, since D(p) is bounded 
from above, nontrivial fluctuations of magnetic held can survive only if D(p) = 0(1 /k) 
for all momenta between the ultraviolet cutoff a -1 and some intermediate scale p, which 
itself is much less than the cutoff. In this case one has < B 2 >oc a _ 4 K _1 and the scale of 
the fluctuations is just right to satisfy eq. ( |3.72| ) in a nontrivial way. 

We see therefore that the emergence of the intermediate scale p is mandatory in 
the compact theory, and the dilute gas approximation for electric vortices gas should be 
reliable. It turns out however, that one can not determine the value of the scale p within 
the dilute approximation itself. The reason is that, at least in the molecular phase, the 
behaviour of the electric vortices is very insensitive to the exact value of the crossover 
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scale /i, as long as it is much smaller than the ultraviolet cutoff. This is so because 
in the Coulomb gas, away from the critical point the vortices are bound in pairs of the 
characteristic size of order of the ultraviolet cutoff a. Therefore if the interaction potential 
is changed at distances x > 1/fi » a, the behaviour of the gas hardly changes at all. 

Even though the precise value of /j, can not be determined, we can establish reliably 
the important fact that /i >> k. This is of course crucial, since if this was not the case the 
propagator would be constant practically for all momenta, and the theory would have no 
propagating degrees of freedom in the continuum limit. The reason this can be established, 
is that, as will be shown in the next section within the dilute gas approximation the scale 
fi is pushed up all the way to the ultraviolet cutoff. Of course, as discussed earlier, one 
cannot believe the approximation for such large [i. However, the dilute gas approximation 
is expected to be valid for values of which are much larger than n as long as they are 
much smaller than the ultraviolet cutoff. 

To verify this picture and also estimate the value of the crossover scale /i we will also 
perform a calculation using a different approximation for the electric vortex gas. This 
approximation is valid for the propagator functions D{jp ) which vanish at large momentum, 
so that functions with the ultraviolet asymptotics of the noncompact propagator can be 
studied reliably. To define this approximation we rewrite the partition function of the 
electric gas in terms of the integer valued field r defined in eq. ( j2.26| ) as 

Z a = j DrDPexp j— TrnK,rp 2 Dr +i J d 2 x2/3(x)(cosr(x) — l)j (3.73) 

Approximating the Lagrange multiplier field (3 by a constant (the same kind of approxi¬ 
mation which lead to eq. ( |3.58|) for magnetic vortices) gives a tractable expression. 

Z a = j Dt exp i^—nnurp 2 Dr +2 a ~ 2 J d 2 x(cosr(x) — l)j (3.74) 
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It is again convenient to define an appropriate fugacity £ in analogy to (|3.69|) . 


n 
47r 


£ = a 2 < cos(t) >= a 2 exp < — 


7T r d 2 p 


n J (27r) 2 - 


Kp 2 D(p) + £ 


-i 


(3.75) 


With the noncompact function D , this fngacity is exponentially small 


£ oc exp(- (an) ) 


(3.76) 


n 

where r is a numerical constant. For this type of variational functions therefore, electric 
vortices can be consistently treated in this “r - approximation”. 


4 Minimization of the energy and the phases of the 
model. 

4.1 The expectation value of the energy. 


Now we are ready to compute the expectation value of the Hamiltonian eqs. fl2.32y2.33|) . 
We start with the magnetic part eq.( [2.32| ). Straightforward integration over the vector 
potential and the noncompact part of the gauge group gives 


< exp{i2ga 2 B(x)} >= exp ^—a i g 2 j ^^p 2 D(p) j x (4.77) 

< exp{i2a 2 7rp(a;)} > p < exp j— 2a 2 n J e tpx D(p)cr(p ) j > a 

Where <> p and <> a mean average over the magnetic and electric vortex ensembles 
correspondingly. The magnetic vortex contribution is trivial, since a 2 p(x ) is an integer. 


We calculate the electric vortex contribution in the dilute gas approximation, keeping 
the first order correction in £. A straightforward calculation gives 

d 2 p 


< exp 2 a n 

[ 87 TK 


exp 


47T 2 

d 2 p 


e ipx D(p)a(p) j > a = 
p 2 D(p) 


(4.78) 


n 


^ + 2C + 2C / d2x [ C0Sh X(x) “ 1 " 2 X(x) 
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where the function X(x) is defined by 


X(x) = —a 2 
n 


d2 P ivx P 2 

4rr2 d^P 2 D ~ 1 (p) + 2 C 


(4.79) 


An analogous calculation for the electric energy part, eq. ( |2.33|) gives 

< exp {^(Xi(x)Ei(x)} >= ex P + ^ 2 D(p)) j (4.80) 

< exp j-^ J ^e ipx p(p)p- 2 D~ l (p)e(p) j > p < exp j ~J ^e zpx a(p)D{p)uj(p) j >* 


where we have defined 


e 


= 7<m, 


(4.81) 


The electric vortex contribution is calculated as 


exp 


47wa 2 


n 


<1 


d 2 p D(p ) 



d 2 ;c 


cosh y(x) 



(4.82) 


with 


Y(x) 


2 f d 2 p ipx p 2 u}(p) 

n J 4tt 2 ^kp 2 /}" 1 ^) + 2C 


(4.83) 


The magnetic vortex contribution this time is nontrivial and is given by 


(7T a 2 

exp <- z 

nn 


d 2 p D 1 (p) 

4vr 2 ^ p 2 D(p ) + 2z 


d 2 x 


cosh L(x) — 1 



with 


L(x) 


2 f d 2 p ipx ejp) 
n J 47t 2 — p 2 D(p ) + 2 z 


(4.84) 


(4.85) 
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Here we also give the expressions for the electric vortex contributions calculated in 


the “r” -approximation: 

< exp 


exp < a g 


2o 2 k J ^e tpx D(p)p 2 r{p) j > CT = 
4 2 f d 2 p KD 2 (p)p 4 /' , 2 


(4.86) 


47t 2 np 2 D(p) + £ 27t 


cosl(x) — 1 — -J 2 (x) 


and 


k 

< exp <- 

7T . 


exp 


2,71Q 2 


n 


f ^e tpx r(p)p 2 D(p)u(p ) 

d 2 p n 2 p 2 D 2 (p) 2 n 

47r 2 np 2 D(p ) + £ 47r j 


>^= 




cos J(x) — 1 + 2 ^ 2 ( x ) 


(4.87) 


with /(x) and J(x) defined by 


/( r ) = 4n q2k [ d2 P c ip* P 2 D(P) 




n J An 2 np 2 D(p) + ( 
2 /' d 2 p ir>x np 2 D(p)uj(p) 


(4.88) 

(4.89) 


47t 2 np 2 D{p)-\-^ 

We will now use these expressions to minimize the average energy. We will consider 
the cases n > 8 and n < 8 separately. 


4.2 n > 


As was discussed in the previous section, for n > 8 the magnetic vortex gas is in the 
dipole phase. The vortex fugacity z vanishes, and in the dilute gas approximation they 
do not contribute to the energy. 

The electric vortex gas will also be treated here in the dilute gas approximation. 
As explained in the previous section, the minimization here is performed on the set of 
variational functions D(p) which have the constant value D(p) = 1/2 k for all momenta 



above some scale /i << a _1 . For these variational functions, the electric vortex gas is 
indeed dilute, and for n > 8 is in the molecular phase. The electris fugacity £ also 
vanishes, and therefore the electric vortices also do not contribute to the energy. 

The calculation is further simplified, since as can be checked explicitly with these 


variational functions D(p) the expectation values in eqs.(4.78) and (4.81) for large n are 


close to one. Therefore the relevant exponentials can be expanded up to first order in 
Taylor series. 

In this approximation, therefore the expectation value of energy is the same as for the 
noncompact theory and the minimization with respect to D(p) gives 



(4.90) 


We therefore find that for n > 8 both magnetic and electric vortices are irrelevant, and 
the compact theory at small momenta is indistinguishible from the noncompact one. 

One assumption that was made in the previous discussion is that the value of the 
propagator D at large momenta is exactly 1/2k. The general argument given in the 
previous section pretty much establishes that this value should be of order 1 /k, but does 
not establish the coefficient. This magic number 1/2 enabled us to conclude that the 
critical value of n for the electric vortex gas is the same as for magnetic vortex gas, n — 8. 
To verify this picture further, we would like to demonstrate explicitly that above the scale 
/i this is indeed the correct behaviour of the propagator. To this end we will minimize 
the energy calculated in the “r-approximation”. This calculation is valid as long as the 
fugacity £ is small. To make sure this condition is satisfied, we will first consider only such 
functions D(p) that behave D(p) oc 1 /\fp z for momenta A 2 < p 2 < a~ 2 with A < a -1 . 
Note, that although A must be smaller than the ultraviolet cutoff, it does not have to be 
parametrically smaller, so that for example A a may remain a finite small constant in the 
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limit a —> 0. 


The terms depending on / and J (see ([4.89Q , (|4.85|) ) in the expression for the expecta¬ 
tion value of the energy can be neglected at large n, since they give corrections supressed 
by powers of 1/n 2 . Now, substituting the formulae of the previous section into equations 
(|2.32|) and (|2.33|) we obtain for the expectation value of the energy H 


2V~ 1 <H>=- 


d 2 p 


D 1 (p) + (p 2 + 4 k )D(p) 


•e 


np 2 D(p ) + £ 


(4.91) 


Minimizing this with respect to D(p) we obtain 

zr 2 G») = z( P 2 + m 2 ) 


e 


( np 2 D(p ) + £) 2 


(4.92) 


where the constants Z and M are 


k <91n£ (np 2 D + ^) 2 f d 2 p 
Z = 1 + — 


^ d In D np 2 D 

M' z = Ak z Z 




K,p 2 D 2 (p) 


np 2 D(p) + ^) 2 


(4.93) 

(4.94) 


The solution of eq.( |4.92| ) is 


D- 1 (p) = Z 1 / 2 (p 2 + M 2 ) 1 / 2 -^p 2 


(4.95) 


Since D 1 is bounded from below, this solution is valid for momenta for which D 1 > 2k, 


that is for 


9 , \t 2 , M 2 . 

p < " 2 = wZ - T> 


(4.96) 


For momenta greater than fi and smaller than A 


D (p > n) — 2 k 


(4.97) 


We see that the physics corresponding to this solution is precisely the same as we were 
describing earlier. Even though we have not determined the value of the scale p yet, let 
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us for the moment assume that it satisfies the relation M « p << a -1 . Then at low 
momenta p 2 « p 2 the variational propagator is the same as in the noncompact case. 
Therefore the physics it describes at all physical scales is the same as in the noncompact 
theory. The only difference is that there appears a nontrivial mass and wave function 
renormalization given by the same factor Z in eqs.( 4.93| , 4.94 ), which is small as long as 
£ is small enough. At scales p 2 > p the propagator is frozen, and therefore at these 
momenta the theory does not contain propagating degrees of freedom. Thus effectively, 
the ultraviolet cutoff has been also changed from a^ 1 to a lower value p. Note also that 
for £ << A 2 , which is where the approximation is valid, the scale p also satisfies p « A. 
The emergence of the scale p and the change in the behaviour of the propagator at large 
momenta is therefore seen already in the regime where the r-approximation is valid. It 
is again difficult to establish the precise value of this scale p. The point is that it grows 
very fast when the auxilliary scale A approaches the ultraviolet cutoff. It is clear from 


eqs. (|4.95| , [4.96|) that the dynamics favours precisely this situation. However, when A is 
close to a -1 , the r-approximation ceases to be valid. 

A rough estimate of p can be obtained by pursuing the r-approximation to the end. 
Calculating the integral in ( |3.75|) using an explicit expression for D(p) ( |4.95|) one has after 
some algebra 

A'k / i y i m \ 

(4.98) 


, 4vr _ 2 1 £ 1 , M 

( = - a exp IIn — 


Using the fact that ln(l/Ma) >> 1 the expression for fugacity is as follows 

47r 


£ = 3 nM In 


na 2 M 2 


[1 + 0 (In ln(l /Ma)/ ln(l /Ma))] 


(4.99) 


From this we get 


p = QnM In 


47T 


na 2 M 2 


(4.100) 
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We see indeed, that /i /M 


o oo. 


We stress again, that the exact value of /j as given in eq. (|4.100|) is not to be trusted 
literally. It is obtained in the region of variational functions D(p) where the nonlinearities 
in the effective r theory are not at all small. For example, there is a large (logarithmic 
in cutoff) renormalization of the fugacity £ due to the cosine interaction already at the 
two loop level. On the other hand we believe that this renormalization will be the leading 
effect. If it is taken into account properly, the expressions ( 4.95| , 4.99| , 4.100) will not be 
changed drastically if £ appearing in them is understood as the complete renormalized 
fugacity. We have checked for example, that the / and J dependent terms in the energy 
expectation value, which appear due to the nonlinearity in the r action, are supressed 
relative to the terms we have kept in eq. ([4.91|) even for fugacities £ of eq. (|4.99|) . We believe 
therefore that the estimate for the crossover scale g eq. (|4.100|) is qualitatively correct. 


To summarize this part, we find that for n > 8 the compactness of the theory is not 
important for infrared physics. The best variational propagator at physical momenta 
is identical to the propagator in the noncompact theory. The only effect of the finite 
radius of compactness, 1/g is that it “freezes” the propagation of the modes with high 
momentum, thereby effectively just changing the ultraviolet cutoff. 


4.3 n < 8 . 


Let us now turn to a more interesting case n < 8. At these values of n both, magnetic 
and electric vortices are in the plasma phase. The fugacities z and £ do not vanish, and 
in the dilute gas approximation we obtain for the expectation value of the energy 
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2V- 1 <H> 


1 

2 


cPp D-'(p) 


(2vr) 2 
+ (p 2 + Ak 2 )D(p) 


2z 


^pW(p) + 2z 


2C 


dr K P 2 D-K P ) + 2C. 


— -^—a 4 C f d 2 x(cosh X — 1 — ^A" 2 ) 

7T J 2 

— — —a~ 2 z [ d 2 x(coshL — 1 — - L 2 ) 

7 T J 2 


(4.101) 


The last three terms appear due to nonlinearities in the 0 and \ Lagrangians. They 
are neglible as long as the appropriate gases are dilute, that is ( and z are smaller than the 
ultraviolet cutoff. We present these terms here for completeness, but will neglect them in 
the following derivations. It can be checked that on the solution we obtain, these terms 
are indeed small. 

Minimizing the energy with respect to D(p) we obtain after some algebra the following 
equation 


2k \ 2 

— p 2 D(p) + 2z 
Tin 

( 1 


27 xnn 


p ip +4 k )+ A( 


V 2nriK 


\ 2 9 

'2k ,, 

J Ak 2 

— p 2 + Bz 
7r n 


(4.102) 


with constants 


A 


1 

167r 4 ?m 


d 2 q _1_' 

47f2 (^■ 1 (?) + 2 0 2 . 


di 2 q q 2 (q 2 + 4:K 2 ) 

^Hdr^D- 1 iq) + 2C) 2 


(4.103) 


and 


B 


4 k 1 f d 2 q 1 

8^n [ ~ Z J A^{^ n q 2 Diq) + 2z) 2 


d 2 q q 2 

47T 2 (^q 2 D (q) + 2z) 2 


(4.104) 
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The solution of this quadratic equation is 

~ 4 ^ ~ <* - 4^)V^)' 


( 2 F(p) 


K 2 p 4 

7 T 2 n 2 


-1 


where we have defined 


, . 4k 2 [% 2 + ^] 

^(p) = i Jj — 1 


(4.105) 


(4.106) 


2 nn^P 2 (P 2 + 4 « 2 ) + ^ 

Let us analyze these expressions. First, note that both values of fugacities are of the 
order of the fugacity of the Coulomb gas. This is true, since as we already know from 
the previous analysis for p 2 > //, the variational propagator is a constant 1/2k. Then the 
integrals in the definitions of the chemical potentials p p and p a in ( |3.65|) and (|3.71|) get 
contributions only from momenta lower than p « a -1 . 

The fugacity of the Coulomb gas behaves like 


— 2 -< 


a e 


(4.107) 


Therefore for n ^ 8 it is smaller, but of the order of the ultraviolet cutoff a 2 . We can 
then estimate the constants A and B as 


A oc 


.-2 


K 


B OC K 


This gives 


4 

K Z 


m oc — v 

a % 


Substituting this into the expression in eq. ([4.105|) we find 


D( 0) = 


0 [za 1 1 1 

- / OC a/t->> - 

£yjF(0) V C « « 2 « 


(4.108) 


(4.109) 


(4.110) 


ft is also easy to see that the function D(p) as given in eq. (14.1051) is a growing function at 
small momenta, and is of the same order of magnitude up to momenta of order p 2 ~ acq -1 . 


34 


















On the other hand the variational function D(p) by definition, can not exceed the value 
1/2k. It follows therefore that the solution eq. (|4.105 ) is unphysical and we have to 
choose for D(p) the end point value 


1 

2k 

at all momenta. 



(4.111) 


We therefore discover that for n < 8 the propagator is constant for all momenta, or 
is completely local in the coordinate space. This should remind one of the behaviour of 
correlation functions in topological theories: there correlation functions of any local oper¬ 
ator are completely local. The situation therefore is such that for n < 8 the Maxwell term 
in our theory becomes totally irrelevant, and the theory degenerates into a pure Chern - 
Simons theory. This conclusion is further confirmed by inspection of the expectation value 
of the energy (|4.102|) . Note, that the structure of this equation is such that for momenta 
which are smaller than z and £ the contributions of the electric and magnetic vortices 
cancell completclly the “noncompact” contributions. Since on our solution Q and z are 
both of the order of the ultraviolet cutoflQ, physical momenta do not contribute at all to 
the energy. This is again consistent with a topological theory, in which the Hamiltonian 
vanishes. 


4.4 Spontaneous breaking of magnetic flux at n < 8. 

Let us discuss in more detail the physical properties of the theory for n < 8. As we have 
just argued, it has many similarities with a pure topological theory. It is clear, however, 
that it can not be equivalent to a simple noncompact pure Chern - Simons theory. First, 

as is obvious from the previous calculation, compactness of the gauge group is crucial 
^Note that with this function D(p) z and ( are both equal to the fugacity of the Coulomb gas. 


35 









in this phase. It is responsible for the appearance of the magnetic and electric vortex 
gases, which are in the plasma phase and thereby determine the infrared properties of the 
model. Second, our derivation in the previous subsection is in large measure independent 
of the assumption k « a _1 . For large k = 0(a -1 ), the “photon mass” is of the order of 
the ultraviolet cutoff, and this case corresponds to the pure Chern - Simons limit of the 
TMGT. It is clear therefore, that the compact pure Chern - Simons theory also undergoes 
phase transition at n = 8, the small n phase being qualitatively different from the large 
n phase. 

We conclude that the compact TMGT at n < 8 is equivalent to compact pure Chern 
- Simons model, but they both must be different from the noncompact Chern - Simons 
theory. What distinguishes the two phases of the pure Chern - Simons theory? 

The fundamental property of a noncompact Chern - Simons theory is the Bohm - 
Aharonov interaction between charged particles. It is interesting to check, therefore 
whether this interaction is still there in the compact theory for n < 8. A straightfor¬ 
ward way to study this question is to calculate the expectation value of the Wilson loop 
in a state which contains a unit external charged. The Bohm - Aharonov interaction 
should show up as the Bohm - Aharonov phase in this expectation value. Introduction of 
a unit external charge at the point x = 0 leads to the following modification of the Gauss’ 
constraint equation 

diEi(x) - 2 kB(x) = gS 2 (x) (4.112) 

The ground state waive functional in this sector should be well approximated by a pro¬ 
jected Gaussian with a nonzero shift in the vector potential. For simplicity, we will take 
the width of the Gaussian to be the same as in the zero charge sector G -1 (p) = k. 

|1 >= J DsiD(f)ex p ji/c [(e„(s* - d'i<j>)Aj + e^d^Sj] - ig<f>(0 ) - ^{Ap 3 - a*) 2 j (4.113) 
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The function a t should be treated as a variational function. Since we are dealing with the 
pure Chern - Simons theory, the function a t should be holomorphic 

a i + id2 = 0 (4.114) 

Additional constraint on a- L follows from the requirement that the state 11 > be normal¬ 
izable. This forces a, to satisfy at zero momentum 

Re eijdidj + Im <9,;a,; = — —6 2 {x) (4.115) 

K 

We will see that the result does not depend on the detailed form of a t . The expectation 
value of the Wilson loop in this state can be calculated in a straightforward manner. 

< exp (ig J BdS) >= exp {ig J eijdidjdS ) < exp J d 2 (j)dS S j ></,< exp n J adS^j > a 

(4.116) 

where for holomorphic functions a.;, the weight for the o averaging is the same as in the 
vacuum state, while the weight for the (j) averaging is given by 

exp{—S 1 ^]} = exp |-^(ai + d^) 2 + ^(a* + a* + <%0 + ie^d^) 2 + ig<p( 0)| (4.117) 

Finally we obtain 

< exp (ig J BdS ^ >= R(S) exp i —^ (4.118) 

where R(S) is a real number. The value of R(S) depends on the phase of the theory, but 
the value of the phase does not. Also, R(S) is nonzero in both phases. 

We conclude therefore that the Bohm - Aharonov interaction of external charges is 
the same in the compact theory for n > 8 and n < 8. 

Nevertheless, the two phases are distinguishible. The operator that distinguishes them 
is not the Wilson loop, but is closely related to its dual. It is the 2+1 dimensional analog 
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of the t’Hooft loop [Q - the operator that creates magnetic vortices. The significance of 
this operator is, that it is an order parameter for spontaneous breaking of the magnetic 


flux symmetry 


10 


In fact, the fate of the magnetic flux symmetry in the compact TMGT is an interesting 
question. Recall, that in noncompact electrodynamics in 2+1 dimensions (both with and 
without the Chern - Simons term) the homogeneous Maxwell equation 


df,Fn = 0 


(4.119) 


ensures the existance of the conserved current. The global charge associated with this 
current is the magnetic flux through the plain 


4 > 



(4.120) 


In QED without Chern - Simons term, this charge is spontaneously broken in the Coulomb 
phase. This breaking is accompanied by the appearance of the massless Goldstone boson 
- the photon. In noncompact TMGT the photon is massive. Accordingly the magnetic 
flux is not broken, but annihilates the vacuum state. In a compact theory the magnetic 
flux is not conserved anymore [|TU|. The magnetic “monopole - instantons” change the 
magnetic flux through the plain by an integer multiple of / g. Consequently, only the 
following subgroup of the flux group remains the symmetry of the theory 


U N = e gN * 


(4.121) 


for integer N. In QED without Chern - Simons term, the magnetic flux is not quantized, 
and therefore the operators Un constitute the group of integers Z. This symmetry again 
is spontaneously broken. Since the group is discreet, its spontaneous breakdown does not 
require an existence of a massless particle, and the photon in compact QED is massive. In 
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TMGT the situation is slightly different. Recall, that the compactihcation of the gauge 
group in TMGT requires the magnetic flux $ to be quantized in units of 2ir/ng. On 
these states only the operators Un for N = 1, ...,n — 1 are represented nontrivially. The 
magnetic flux symmetry in compact TMGT is therefore just Z n . Let us note, that i t is 
this reduction of the magnetic flux group to Z n that is responsible for the mixing of the 
states with n/2 positively and negatively charged particles in the spontaneously broken 
SU (2) TMGT discussed in paper |lT| . 

The natural question is, what is the realization of this Z n flux symmetry in compact 
TMGT. To answer this question we may calculate the vacuum expectation value of the 
appropriate order parameter. This order parameter should be an operator which is gauge 
invariant under both, compact and noncompact gauge groups, and should transform non¬ 
trivially under the flux symmetry. The suitable operator is 

. 27 rm 


v m (x, C ) = exp l i f 
an Jc 


eijdxiEj 


(4.122) 

gn Jc " 

The contour C is a semiinfinite line with an endpoint at x. For integer m < n this 
operator commutes with the elements of the compact gauge group. The calculation of 
the expectation value of v is straightforward and parallels exactly the calculation of the 
electric part of the expectation value of the energy. Without giving the details here, we 
just describe the result. 


For n > 8, there are no vortex contributions and the result is basically the same as in 
the noncompact theory. It has the form 

< v m >= exp{— Kar l L} (4.123) 

where K is a numerical constant, and L is the length of the curve C. In the infinite 
volume, L —> oo, and we find 

< v >= 0, n > 8 (4.124) 
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For n < 8 the contributions of the magnetic and electric vortices cancell exactly the 
noncompact contributions, precisely in the same way as they did in the calculation of the 
energy expectation value. There is therefore no linear divergence in the infinite volume 
limit, and we obtain 

< v 0, n < 8 (4.125) 

We conclude therefore, that for n > 8, the Z n magnetic flux group in compact TMGT 
is unbroken, and for n < 8 it is broken spontaneously. These different ways the symmetry 
is realized in the vacuum distinguish between the two phases of the model. 

5 Discussion. 

To summarize, the result of our variational calculation is the following. For n > 8 the 
magnetic and electric vortices in the wave function are bound in pairs. Their only effect 
is to introduce an intermediate scale /i, below which the physics is the same as in the 
noncompact theory. This scale // becomes infinite in the limit of inhnte ultraviolet cutoff, 
and therefore in the continuum limit the effects of compactness disappear. This result is 
in qualitative agreement with the results of previous studies ||. Here we want to make 
the following remark. The basic picture of ref. |3] is that the monopoles in the compact 
Chern - Simons theory are bound in pairs by a linear potential. On the other hand 
the magnetic (and electric) vortices in our ground state wave functional have interaction 
which is logarithmic at large distances. It is important to realize that these two claims 
are not inconsistent with each other. The point is, that the object which in path integral 
formalizm is represented by a monopole - antimonopole pair has a segment of a Dirac 
string, stretched between them. This Dirac string becomes “observable” in the Chern - 
Simons theory and carries a finite action density, which is the origin of the linear potential 
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between the monopoles. In our Hamiltonian description, a widely separated monopole 
- antimonopole pair is represented by a configuration in the wave function which has a 
nonvanishing magnetic vortex density p at two points, but also a string of electric vortex 
dipoles cr along the line that connects these points. This line of electric vortex dipoles 
corresponds precisely to the “observable” segment of the Dirac string. Configurations of 
this type are indeed suppressed in our wave functional by an exponential of the length 
of the dipole string segment. In this sence it is indeed true, that monopoles in our wave 
functional interact with a linear potential. 

The new result we find, is that for n < 8 the nature of the ground state of the theory 
is very different. Due to the liberation of magnetic and electric vortices, the correlation 
functions of local observables (such as B or E) become completely local. The theory 
therefore does not describe any propagating degrees of freedom. The liberation of the 
vortices also leads to spontaneous breaking of the magnetic flux symmetry in this phase. 
In this respect this phase is similar to the confining compact QED without the Chern - 
Simons term, whose vacuum also spontaneously breaks magnetic flux. 

The dichotomy between the “monopoles” on one hand and the Chern-Simons term on 
the other hand, is resolved therefore in this nontrivial way. 

Finally, let us remark, that in recent years TMGT have seen many applications in 
condensed matter physics, both in relation to Quantum Hall Effect and high temperature 
superconductivity. In the latter case a prominent role is played by the so-called semions. 
Those are charged particles coupled to Chern-Simons theory with n = 4. It is clear from 
our results, that the question of compactness of the Chern-Simons theory must be crucial 
for the “semion physics”. The physics usually discussed, corresponds to the noncompact 
theory. In the compact case the semions must have a very different behaviour. We hope 
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to return to this question in future work. 
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